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(57) Abstract 

me invention relates to a method and a 
is no. based on the assumption of prior art methods that the P^Kdnft * ukw ™" modu)ator ^ resolulion f | e ments of 

imaging operation with a simple one-to^e , mapping between resdunon element ota P fln jmcnsity patlern wilh , w 

the ge^rated intensity panem is provded. ^^^ ,B ^^nS rXtion w.th a spatial phase mask for modulat.on of the 
loss of electromagnetic energy, compnsmg spat.al on Sa ^solution e.ements of the spatial phase mask, each phasor 

phase of the incident eletromagnetic rad.at.on by « . atujn prede , e nnined values for predetermined 

value being determined in such a way that me values o^JSJ ™ s SfpaSphase mask corresponds to a dist.nc. .ntcnsit, level 
spat.al fluencies, and the phasor value of the specific ^^7"'° ^nner for phase shifting of a pan of the cleetromagntic 
oV the .mage of the resolution element m the mtens.ty P™™-"* '^IVv panem bs m.erference .n the .mage plane of the .mag.ng 
rad.at.on. m comb.nat.on *.th an .magmg system for gcnenmono JJJJ ^ „„„ ^ thf wmaini „ F pan of .he 

ssstem between the pan of the electromagnet* ra.dat.on that has been phase srm.eu r 
electromagnetic radiation 
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Phase Contrast Imaging 
' FIELD OF THE INVENTION 

The invention relates to a method and a system for synthesiz- 
ing a prescribed intensity pattern based on phase contrast 
imaging . 

5 BACKGROUND OF THE INVENTION 

It is well known to form an image on an illuminated surface 
of a body by absorption or blocking of energy of an illumi- 
nating beam. For example in an overhead projector, an over- 
head transparent absorbs or blocks part of the light beam of 
10 the projector whereby a large image of an overhead is formed 
on a screen. However, this results in a loss of light inten- 
sity as part of the emitted light from an image forming 
system is reflected or absorbed. 

To avoid loss of energy causing, e.g. loss of light intensity 
15 of the synthesized intensity pattern, power dissipation 

generating heat in components of the system, etc., methods 
and systems have been developed wherein the phase of a light 
beam is modulated instead of the amplitude or intensity of 
the light beam, as modulation of the phase of the light beam 
20 do not lead to loss of energy. The phase modulation is fol- 
lowed by a conversion of the phase modulation into an ampli- 
tude or intensity modulation. 

A diffractive optical element, such as a holographic optical 
element, may be used to generate a phase modulation. Then, 

25 the resulting intensity modulation at each point of a picture 
formed by conversion of the phase modulation into intensity 
modulation will depend upon the phase modulation values at 
each point of the diffractive optical element as the light 
intensity at each point of the picture is formed by a coher- 

30 ent superposition of light received from the entire surface 
of the diffractive optical element. Diffractive optical 
elements are rather complex to design for synthesis of a 
prescribed intensity pattern. 
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tmaging methods and system may also be used in «—c"- 
„ith phase adulation. These methods and system are charac- 
terised by Che fact that the intensity of a point of a pic- 
ture formed by conversion of phase Elation into ntensity 

; modulation will depend upon the phase modulation value of one 
point of the phase modulator only as this point is imaged 
onto the picture point in question by the imaging system 
This one-to-one relationship maxes the design of phase modu- 
lators in these systems simple. Methods and systems of this 

0 ki nd are named phase contrast imaging methods and systems. 

Phase contrast imaging methods were originally developed 
within the field of microscopy. Many objects «* ^ 
microscopy are largely transparent, thus absorbing little or 
no light. When light passes through such an object, the 
5 predominant effect is the generation of a spatially varying 
Phase shift which can not be seen by a human as the eye of a 
human responds to light intensity and colour and does not 
respond to the phase of light. 

in 1935. Fritz Zemik proposed a phase contrast technique 
20 which rests on spatial-filtering principles and has the 

advantage that the observed intensity is linearly related to 
the phase shift introduced by the object. 

Suppose that a transparent object with amplitude 
transmittance 

25 t(x,y> - explj*(x,y)I 111 

is coherently illuminated in an image- forming 
simplicity, a magnification of unity is assume and th 
finite extent of the exit and entrance pupils of the system 
is neglected. Further, a necessary cond "ion to achieve 
30 linearity between phase shift and intensity is that h phase 
shift * be less than 1 radian, in which case the amplitude 
transmittance can be approximated by 
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t(x,y) = 1 + j<t>(x,y) 



(2) 



The terms of order 4> 2 and higher are neglected in this appro- 
ximation. It is seen that the first term of (2) leads to a 
strong wave component that passes through the sample without 
5 change, while the second term generates weaker diffracted 
light that is deflected away from the axis of the system. 

The image produced by a conventional microscope can be writ- 
ten 

I » | 1 + j*l 2 - l (3) 

10 where the term <*>* has been approximated by zero. It is seen 
that the diffracted light is not observable because it is in 
phase quadrature with the strong background. As Zernik recog- 
nized that the background is brought to a focus on- axis in 
the focal plane while the diffracted light - containing 

15 higher spatial frequencies - is spread away from the focal 
point, he proposed that a phase -changing plate be inserted in 
the focal plane to modify the phase relation between focused 
and diffracted light. 

The phase- changing plate can consist of a glass substrate on 
20 which a small transparent dielectric dot has been coated. The 
dot is placed at the center of the focal plane and has a 
thickness and index of refraction such that it retards the 
phase of the focused light by either tt/2 radians or 3ir/2 
radians relative to the phase retardation of the diffracted 
25 light. In the former case the intensity in the image plane 
becomes 

I = |exptj(ir/2) + j4>| 2 - |jd+*>l 2 - 1 + 2 * (4) 
while in the latter case 



I - |exp[j(3*/2)+j<*>| 2 = |j(-W)| 2 -1-2* 



(5) 
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Th us, the image intensity has become linearly 

ph ase shift ♦. When the phase of the background -^"^ 

/, P h. result is known as positive phase contrast, wnn 
IVs/^etaraatlon is said to yield negative phase contrast. 

• it is seen that the method described above leads to a phase 
' n rast Uing method that provides a small phase srgnal 
that is superimposed on a large DC-component. This leads to 

it will be necessary to attenuate tne 

„. „ contained in the phase modulated signal. 
0 the information contained in cne p 

However, the attenuation of the ^ 0m ^no^ ^ Field 
energy. This kind of filtering is usually denoted Dark 



Filtering . 



15 
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it is another disadvantage of the phase contrast imaging 

described above that it is based on the assumption 
^tlht ^a e shift * is less than 1 radian which is very 
that the P practical real-life applications, 

often not » P such applicat ions. 

However, the " ^J^ssic action is not ful- 

rrrr, sirr^i— technical 

,» EP 0 657 760 a phase contrast imaging system is 

^ 11 an image simulation and ^-^^^ 

on the Texas Instrument flexure beam digrta mirror dev 

-L- ^ refLrd ^r«r,r^r- is 

ilging method. The flexure beam DMD provides a nicker 
Elated wave and accordingly ..tica 

, ssr^ r ■ incxudes che 

corresponding disadvantages described above. 

of a Phase contrast imaging system is 
T dTTcB 2 ». 7i6. -herein an optica! guide-beam 

£££ - " — syscem is disclosed that 
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provides a spatially intensity modulated guide-beam. A spa- 
tial phase modulator is used to generate the guide-beam. The 
phase encoding of the spatial phase modulator constitutes a 
Periodic square-wave modulation (50% duty cycle) of two phase 
3 values 0 and */2. The phase modulation is converted into an 
amplitude modulation by Fourier transforming lenses and a 
phase plate providing a phase shift of the background signal 
by tt/2. A method for synthesizing the specific intensxty 
pattern of the optical guide-beam based on phase contrast 
imaging is not disclosed in this document. 

A similar example of a phase contrast imaging system is 
disclosed in "Array illuminator based on phase contrast", 
Applied Optics Vol. 27, NO. 14, pp. 2915-2921 (1988). A 
method is disclosed of converting a wide beam of uniform 
intensity into an array of bright spots without losses. The 
input spatial phase mask constitutes a periodic array of 
phase dots with the phase value *, the remaining area of the 
phase mask having the phase value 0. The phase modulate xs 
converted into an amplitude modulation by Fourier 
transforming lenses and a phase plate providing a phase shxft 
of the background signal by *. The method is limited to the 
implementation of periodic array configurations with the 
binary phase values 0 and it. 

It is well-known to nse so-called "radiation f ocusators" . 
25 i.e. counter generated holographic optical elects for 

spatial phase modulation of a light beam, e.g as drsclosed in 
special issue on Computer Optics in the USSR. Optics and 
Lasers in Engineering. Vol. 15. no. 5 1991. However, such 
elements are complicated to synthesize. Typically, they are 
30 synthesized in such a way that the desired image is formed xn 
the Fresnel region or the Frauenhofer region. Thus, the 
intensity of a resolution element in the generated image 1S a 
function of several, typically all, phase values of the 
resolution elements of the holographic optical element. 
35 Obviously, this complicates the design of a general purpose 
holographic optica! eiement and advanced, very time consume 



20 
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i« M have to be applied. Further, the complicated 

ISlt Sii *— "V -ngeable 3P«-1 phase 
modulators with such elements. 

6 Tt is a further disadvantage of holographic optical elements 
5 L a carr^r frequency is needed to 

lig ht from non-diffract- light - « _ 

system geometry and a need for a arrir 

support these high frequency terms. 

10 SUMMARY OF THE INVENTION 

It is an object of the present invention to provide an appar- 
atus of the above *ina which apparatus is robust, compact, 
Zle to design and relatively cheap to manufacture. 

nf the oresent invention to provide an 
15 ZZZZX aPpTtus for phase contrast i^ging that 
take all terms of the Taylor's series: 

v <b{x,y) 2 ,*U 1 yli + ll£EiZ^ + . .. (6) 
tix.y) =e^ (x ' y) -i+i*U.y» ' 2! — 3! ^ 

-v,,,* is not based on the assumption that 
into account and^ ° ^ 1 radian. It is i^ortant to 

20 n of th. ; Taylor, series contribute 

~ „ be. «, * 

represented by numeral 1 in equation (2) . 

It is still another object of the present invention tc .pro- 
's - - -rrrre d ed a :: rr^zr c£ 

r/Snrr—eThe infection contained in the phase 
modulated signal. 
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It is yet another object of the present invention to provide 
an improved method based on a simple imaging operation with a 
simple one-to-one mapping between resolution elements of a 
spatial phase modulator and resolution elements of the gener- 
5 ated intensity pattern. 

According to the invention a method is provided for syn- 
thesizing an intensity pattern with low loss of electromag- 
netic energy, comprising spatial modulation of electromag- 
netic radiation with a spatial phase mask for modulation of 
10 the phase of the incident electromagnetic radiation by phasor 
values of individual resolution elements of the spatial phase 
mask, each phasor value being determined in such a way that 

1) the values of the Fourier transformed phasors attains 
predetermined values for predetermined spatial fre- 
15 quencies, and 



20 



2) 



the phasor value of a specific resolution element of 
the spatial phase mask corresponds to a distinct 
intensity level of the image of the resolution 
element in the intensity pattern, 



and a spatial phase filter for phase shifting of a part of 
the electromagnetic radiation, in combination with an imaging 
system for generation of the intensity pattern by interfer- 
ence in the image plane of the imaging system between the 
part of the electromagnetic radiation that has been phase 
25 shifted by the phase filter and the remaining part of the 
electromagnetic radiation. 

Although, the present method is related to encoding of spa- 
tial phase masks in two spatial dimensions (planar encoding) , 
the principles of the method may be utilized for phase encod- 
30 ing in one to three spatial dimensions and/or in the temporal 
dimension. 
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The electromagnetic radiation may be of any frequency range 
The electron^ frequency 
of the electromagnetic spectrum, i.e. the gamma h 

n ae the ultraviolet range, the visible range, the infrared 
range, the ultra range, the microwave 

i-v,o far infrared range, the X-ray range, 
range, the far 3 oresent method 

tip (hi ah frequency) range, etc. The present 

TZ'soZZ^ » — — ■ such 38 electron 

radiation, neutron radiation, etc. 
P^ly. - -tro^t- — 

"e tn tensity pattern is synthesis by —e ° » 
electromagnetic waves emitted from a co-=» source of elec 
trlmagnetlc radiation but the phases of which have been 
^differently, it is retired that the f reouency range 

7"::::* ZZZZ ^their superposition aerates the 

T f the freouency range is too 
desired intensity pattern. If the freque y 

nroad the two waves will be incoherent and the phase mfor 

. in :» be iost as r:. 

: etromagnetic radiation to be superpos itioned is less 
th an the wavelength of the radiation. Thi. , l * 

i-hat allows the electromagnetic radiation to u 

l TZ Z a Hg-lamp can be used as a light source in a 
s racco^g tTthe present invention ----- 
h,» s laser light source that the speckle noise is 
"Teed L requirements of the spatial coherence of the 
30 :rtrl^rradiation depend upon the space bandwith 

Luct of the corresponding system -J^^T^. 
required system performance is to the theoretic 
able performance of the system. 
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• Preferably, the electromagnetic radiation is generated by a 
coherent source of electromagnetic radiation, such as a 
laser, a maser, a phase-locked laser diode array, etc. How-^ 
ever a high pressure arc lamp, such as a Hg lamp, a Xe lamp, 

5 etc may also be used and even an incandescent lamp may be 
used as a source of electromagnetic radiation in a low per- 
formance system. 

A spatial phase mask is a component that changes the phase of 
an electromagnetic wave incident upon it. The spatial phase 
10 mask may transmit or reflect the incident electromagnetic 
wave. Typically, the spatial phase mask is divided into a 
number of resolution elements each of which modulates the 
incident electromagnetic wave by changing its phase by a 
specific predetermined value. The predetermined values are 
IS assigned to each resolution element in different ways depend- 
ing upon the technology applied in the component. For example 
in spatial light modulators, each resolution element may be 
addressed either optically or electrically. The electrical 
addressing technique resembles the addressing technique of 
20 solid-state memories in that each resolution element can be 
addressed through electronic circuitry to receive a control 
signal corresponding to the phase change to be generated by 
the addressed resolution element. The optical addressing 
technique addresses each resolution element by pointing a 
25 light beam on it, the intensity of the light beam ««™PO-- 
ding to the phase change to be generated by the resolution 
element illuminated by the light beam. 

Spatial phase masks may be realized utilising fixed phase 
»SKS. devices comprising liquid crystals and being based on 
30 liquid crystal display technology, dynamic mirror devices, 
digital micromirror arrays, deforce mirror devices, mem- 
brane spatial light modulators, laser diode arrays (inte- 
grated light source and phase modulator), smart pixel arrays. 



etc, 
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A spatial phase filter is typically a fixed phase mask, such 
LTn o tically flat glass plate coated with a - - ~ 

iayer at specific — ^ sec^on J« 

spatial phase masks mentioned xn the previous 
5 also be used for spatial phase filters. 

The imaging system maps the phase modulating resolution 
laments of the spatial phase mas, on the targe, surface of 
! h e synthesized intensity pattern. It may comprise a 4f 

conflation .two Fourier transforming lenses uti izmg 
10 lesion of light or one Fourier transforming lens u « 

izing reflection of light, or a single imaging lens. However. 

any optical imaging system providing a filtering plane for 

Z spatial phase filter may he applied in a phase contrast 

imaging system. 

15 in the method according to the present invention, the syn - 
Resized intensity pattern is generated by »*«»~"~ * 
Z electromagnetic waves in the image plane of ^ imaging 
Z r *n The spatial phase mask changes the phase values of an 
system. The spatial v .„„„,, fhe imaaing system 

electromagnetic wave incident upon it and the imaging 
20 directs the electromagnetic wave with changed phases 
20 r e lected from or transmitted through the spa, :ial phase mask 
cowards the spatial phase filter. The phase filter phase 
l*«ts a part of the electromagnetic radiation and the imag 
n systems adapted to superimpose in the isage plane the 
25 Sse shifted part of the electromagnetic radiation with the 
part of the electromagnetic radiation that is not phase 
shifted by the spatial phase filter. 

w-,. rte soatial Phase filter is positioned in the 
30 a lens while the spatial p ele ctromagnetic 
back focal plane of the lens, wnereoy 

field at the phase mask is Fourier transformed by th lens 
into a second electromagnetic field at the phase filter. 
Thu . s cific spatial freguencies of the first electromag. 
35 Lie field will be transmitted through the spatia! phase 
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filter at specific positions of the phase filter. For 
instance, the energy of the electromagnetic radiation at zero 
regency (DC) is transited through the phase filter at the 
injecting point of the Fourier plane and the 
5 of the lens also denoted the zero-order diffraction region. 

lt is presently preferred that the spatial phase filter is 
adapted to phase shift the DC-part of the electromagnets 
Edition and to leave the regaining part of the 
netic radiation unchanged or. alternatively, to leave the DC 
L0 part of the electromagnetic radiation unchanged and to phase 
shift the remaining part of the electromagnetic radiation 
The last alternative is preferred when the energy level of 
the DC-part of the electromagnetic radiation is so high that 
the phase shifting part of the phase filter will he destroyed 
15 by it. For example in laser cutting, the DC level of the 
laser beam can be so high that a phase shifting dot posi- 
tioned at the intersecting point of the DC part of the laser 
beam at the phase filter would evaporate. It is ^« 
co block the electromagnetic radiation (no transmittance) in 
20 the zero-order diffraction region, however, the DC energy of 
the radiation is then lost. 

Below, an expression of the intensity of the synthesized 
intensity pattern as a function of the phasor values »(>,y) 
of the phase mask, when the DC-part of the electromagnetic 
25 radiation is phase shifted, is deduced. 

Electromagnetic radiation incident on the spatial phase mask 
can be described by a function A(x,y) . where A(x.y) is a 
convex number (amplitude and phase, of the incident field on 
the point (x,y» of the spatia! phase mask. At the point 

30 <x,y> . the spatial phase mask modulates the phase of the 
incident radiation with a value *,x,y, so that the 
after reflection by or transmission through the spatial phase 
mask may be described by the function A(x.y) * e • . 
e i*<*.r> being the phasor value of the point (x.y) of the 

35 spatial phase m*sk. As A.x.yl preferably is a constant vaiue 
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• over the entire surface of the spatial phase mask, the term 
is left out of the following equations for simplicity. 

The expression of the electromagnetic radiation incident on 
the spatial phase filter may now be separated into an AC- term 

5 and a DC- term. If the DC- term of the field is denoted a. the 
AC-term of the field is given by the term e lo(x ' y) -a. As the 
spatial phase filter changes the phase of the DC-part of the 
electromagnetic radiation by d, the intensity of the syn- 
thesized intensity pattern at the image plane of the imaging 

10 system is given by: 



I(x',y') - \e^ (x '' y,> + a(e i0 -D| : 



(7) 



wherein (x',y') is the coordinates of the image of the point 
(x,y) of the spatial phase mask formed by the imaging system 
in the image plane. 

15 It should be noted that the second term of the equation is a 
complex number that adds to the phasors of the spa- 

tial phase mask and may be interpreted as a contrast control 
parameter for the synthesized intensity pattern I(x',y'). 

According to a preferred embodiment of the invention, the 
20 average value of the phasors is adjusted in order to control 
the range of intensity levels. 

instead of phase shifting the DC-part of the electromagnetic 
radiation, it is also possible to synthesize a prescribed 
intensity pattern by phase shifting other parts of the elec- 

25 tromagnetic radiation by adapting the spatial phase fxlter to 
phase shift electromagnetic radiation incident upon one or 
more arbitrary regions of the phase filter and leaving the 
phase of the remaining part of the electromagnetic radiation 
unchanged and then superimposing the two parts of the elec- 

30 tromagnetic radiation. The corresponding mathematics and the 
corresponding design procedures for the spatial phase mask 
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and spatial phase filter will of course be more complicated 
than for the method described in the previous section. 

A simple example of phase shifting a part of the electromag- 
netic radiation of a spatial frequency different from the 

5 zero frequency is provided by moving the DC-part of the 

electromagnetic radiation to another spatial frequency in the 
Fourier plane (identical to the plane of the spatial phase 
filter) utilizing an optical component with an appropriate 
carrier frequency (i.e. a grating or a prism) or, preferably, 

10 encoding the function of a grating or a prism into the spa- 
tial phase mask, and adapting the spatial phase filter to 
change the phase of the electromagnetic radiation at this 
spatial frequency and to leave the phase of the remaining 
part of the electromagnetic radiation unchanged. 

15 According to another preferred embodiment of the invention, 
the phase mask is not positioned in the back focal plane of 
the lens but in the Fresnel region of the lens instead. In 
this case, the electromagnetic field at the phase filter will 
be given by a Fresnel transformation of the electromagnetic 

20 field at the spatial phase mask. This further complicates the 
mathematics and the design procedures, for example the term a 
in equation (7) has to be substituted by the value of the 
Fresnel transformation at the point (s) of phase changes of 
the phase filter. However, the Fresnel transformation may be 

25 calculated from a Fourier transformation by multiplication of 
the phasor values of the spatial phase mask by a quadratic 
phase factor followed by a Fourier transformation. 

It is an important aspect of the present invention that each 
intensity level of the synthesized intensity pattern for each 
30 resolution element may be generated by at least two different 
phasor values of a resolution element of the spatial phase 
mask. 
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For example, when the spatial phase filter phase shifts the DC 
part oTthe electromagnetic radiation, it will he shown later 
that, advantageously, the average « of the phasors of the 
resolution elements of the phase mask should be equal to tf 
; and the value of the phase shift 6 should be equal to *. In 
this case, the intensity of the synthesized image pattern at 
the image <x',y') of the resolution element (x,y) will be 
given by: 



l(x',y') = 2(1 - coe<t>(x' ,y')) 



(8) 



LO It is seen that complex conjugate phasors (values of + of 

opposite sign) result in identical intensity levels I(x ,y ). 
It can be shown that for any value of the modulus of the 
average of the phasors | a| , two phasors exist that will 
generate identical intensity levels of the synthesized inten- 

15 sity pattern. 

Purther. if the spatial phase filter phase shifts part. .of 
the electro^gnetic radiation different from the DC-part the 
phasor value that generates a specific intensity level will 
depend on the position of the resolution element xn question. 
20 i.e. the phasor value and the position of the resolution 

element with that phasor value together define the intensity 
level at the image of the resolution element in the syn- ^ 

thesized intensity pattern. Still, it is true that for each 

resolution element of the spatial phase mask, each intensity 
25 level of the synthesized intensity pattern may be represented 

by one of two different phasors of complementary phase 

values . 

This freedom of being able to select, for each intensity 
level to be generated and for each resolution element of the 
30 spatial phase mask, one of two phasors is used to control the 
phase of the'Fourier transform of the phasors at specie 
spatial frequencies by selection of phasors with appropriate 
phase values to ensure two intervals of biunique functional 
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dependence 
values . 



between phasor values and corresponding intensity 



This freedom of choice of phasors may be utilized to select 
phasors of neighbouring resolution elements of the spatial 

5 phase mask with a maximum difference between them, thereby 
generating an electromagnetic radiation emitted from the 
phase mask with a maximum content of high spatial frequencies 
which will generate a good separation of the DC part of the 
electromagnetic radiation from its AC part. However, any 

10 other strategy of selecting between two possible phasor 

values of each resolution element may be chosen to generate a 
desired spatial frequency content of the electromagnetic 
radiation. 

Preferably, the phase of the Fourier transform of the phasors 
15 at specific spatial frequencies is adjusted in order to 
control whether the relation between each phasor and the 
corresponding intensity level is a monotonic increasing or a 
monotonic decreasing function. 

Below, a set of different methods are described that are 
20 provided according to the present invention for adjustment of 
the modulus of the Fourier transform of the phasors at speci- 
fic spatial frequencies to attain a prescribed value. If 
convenient, the methods may be combined. 

According to one of the methods, the individual phasors of 
25 the resolution elements of the phase mask are adjusted by a 
constant value until the desired value of the modulus of the 
Fourier transform of the phasors at specific spatial fre- 
quencies is attained while maintaining prescribed relative 
intensity levels between intensities of resolution elements 
30 of the intensity pattern, i.e. iteratively. 

According to another method, the individual phasors of the 
resolution elements of the phase mask are adjusted utilizing 
histogram techniques known from image processing. A histogram 
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v, , , hart showing the number of resolution elements of 
is a bar chart showing u sriec ific intensity 

= redetermined distribution, etc., may u= 
i histogram to a predetermi „ ourier transform of the 

iterative!, until the modulus of 'he Fourier tra 
phasors at specific spatial frequencies attain P 



value . 



10 



15 



rrS^tSSL - the phasors at specific spatial 
frequencies. 

modulation, etc. 

r n^sTcLUce of the " 

0£ the Fourier transform «£^cZ£T*>^ of 
frequencies. Thus, it is preferred to of 
the radiation source 1^ inten- 
the intensity pattern so that a se 4 
5 sity patterns show uniform intensity levels. 

-« » oref erred embodiment of the invention, the 
According to a P" f «™ tQ ^ the spatial 

8haP e of the phase i ^ ^ phase 
frequency content of the pnaso f che 

e g. to optimize the desired separation of the part 
,0 electromagnetic radiation to be phase filtered from the 
remaining part of the el ectrom*gnetic radiation. 
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It is within the scope of the present invention that the 
imaging system further comprises zooming means for variable 
scaling of the synthesized intensity pattern. The zooming of 
the imaging system may be dynamically controllable, e.g. in 
5 response to the scaling of the pattern of phasor values of 
the phase mask. 

According to the present invention, the power of the radi- 
ation source may be controllable in response to the spatial 
scaling of the pattern in the phase mask and/or the zooming 
10 of the focusing system. 

In order to provide a compact and integrated system according 
to the present invention, the optical function of a Fourier- 
transforming lens is encoded into the phasors of the spatial 
phase mask. The Fourier transforming lens may be refractively 
15 or dif f ractively encoded into the phase mask. 

Similarly, the optical function of an output lens may be 
encoded into the phase filter either refractively or 
dif f ractively . 

Further, a compensation may be encoded into the phasor values 
20 of the spatial phase mask so that part of the electromagnetic 
radiation modulated by the phase mask has a substantially 
flat intensity profile in the image plane. Without this 
compensation, part of the electromagnetic radiation modulated 
by the phase mask will have a flat profile with perturbations 
25 resulting from the phase filtering superpositioned upon it. 
This may cause "ringings' 1 (oscillations) at the edges of the 
synthesized intensity pattern. 

According to another preferred embodiment of the invention, 
the source of electromagnetic radiation comprises one or more 
30 light sources of different wavelengths corresponding to three 
different colours, such as red, green and blue, for gener- 
ation of intensity patterns of arbitrary colours. Further, 
several independent systems each one illuminated by its own 
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' wavelength can be combined into a single multi-wavelength 
system. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig . i shows a 4f optical system for phase contrast imaging. 

5 Fig. 2 shows a 2f optical system for phase contrast imaging, 

Fig. 3 shows a If optical system for phase contrast imaging, 

Fig. 4 shows (A) off-axis read-out of reflective SLM and (B) 
on-axis read-out of reflective SLM. 

Fig. 5 shows schematically an example of a prescribed inten- 
10 sity pattern in ID. 

Fig. 6 shows schematically the resulting phase encoding 
corresponding to Fig. 5. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Pig i shows a 4f phase contrast imaging system (1) • A laser 
15 «) emits a light beam which is expanded by a beam expander 
3 into a plane light wave of uniform intensity and directs 
It towards a spatial phase masK ,4, . The light beam « trans- 
»itted through the spatial phase mas* («) and a 
transforming lens (5. . The spatial phase masK rs P°«"°^ 
20 in the front focal plane of the lens (5) and a spatial phase 
filter («) is positioned in the back focal plane o« th< ^lens 
,5, that is also the front focal plane of a lens (7) . The 
Pourier transforming lenses .5. 7, need not have identical 
£ ocal lengths. Different focal lengths lead to a magnifies- 
25 tion ratio different from one. The phase filter CO phase 

shifts the zero order diffraction part <8> of the light phase 
adulated by the spatial phase mas* <4, . The synthesize 
intensity pattern is generated in the bacK focal plane , 



(9) of 
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the lens (7) and a dynamic focusing system (10) images the 
synthesized intensity pattern onto a focusing plane (11) . 

The optical system is controlled by a computer (12). The 
computer (12) comprises interface means for addressing each 
5 of the resolution elements of the phase filter (4) and trans- 
mitting a phasor value to the addressed resolution element. 
Further, the computer (12) comprises laser control means for 
controlling the power of the laser (2) and imaging control 
means for controlling the focusing and the image ratio of the 
10 dynamic focusing system (10) . The computer (12) also com- 
prises input means, such as a keyboard, a diskette drive, an 
optical disc drive, a network interface, a modem, etc, for 
receiving an image pattern to be synthesized by the system 
(l) From the received image pattern, the computer is adapted 
15 to calculate phasor values to be transmitted to the resol- 
ution elements of the phase mask, e.g. based on a histogram 
technique as described herein. Optionally, the phase shift of 
the phase filter (6) is adjustable and controllable by 
optional phase control means of the computer (12) which may 
20 be further adapted to adjust the phase shift, e.g. utilizing 
equation (18) . 

Fig. 2 shows a 2f phase contrast imaging system (20) . A laser 

(21) emits a light beam which is expanded by a beam expander 

(22) into a plane light wave of uniform intensity and directs 
25 it towards a spatial phase mask (23) and a polarization beam 

splitter (24) and a quarter-wave plate (25) . The polarization 
beam splitter (24) and the quarter-wave plate (25) allows 
beam- splitting of light of a specific linear polarization 
without the power loss associated with conventional beam- 

30 splitters due to splitting of the beam in both directions of 
transmission through the beam-splitter. After transmission 
through the polarization beam splitter (24) and the quarter- 
wave plate (25), the light beam is transmitted through a 
Fourier transforming lens (26) and is reflected from a spa- 

35 tial phase filter (27). The spatial phase mask (23) is posi- 
tioned in the front focal plane of the lens (26) and the 
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' spatial phase filter (27) is positioned in the back focal 
plane of the lens (26) . The phase filter (27) phase shifts 
the zero order diffraction part (28) of the light that is 
phase encoded by the spatial phase mask (23) . The synthesized 

5 intensity pattern is generated in the back focal plane (29) 
of the lens (26) and a dynamic focusing system (30) images 
the synthesized intensity pattern onto a focusing plane (31) . 
As described for the system shown in Fig. 1. the system (20) 
is controlled by a computer (32) . 

10 Fig. 3 shows a If phase contrast imaging system (40) . A laser 

(41) emits a light beam which is expanded by a beam expander 

(42) into a plane light wave of uniform intensity and directs 
it towards a spatial phase mask (43). The light beam is 
transmitted through the spatial phase mask (43) and an image 

15 forming lens (44) . A phase filter (45) positioned in the back 
focal plane of the lens (44) phase shifts the zero order 
diffraction part of the light phase encoded by the spatial 
phase mask (43) . The synthesized intensity pattern is gener- 
ated in the image plane (46) of the lens (44) and a dynamic 

20 focusing system (47) images the synthesized intensity pattern 
onto a focusing plane (48) . As described for the system shown 
in Fig. 1. the system (40) is controlled by a computer (49). 

Fig. 4 shows details of (A) an off-axis read-out of a reflec- 
tive phase mask (50) (or a spatial light modulator) and of 
25 (B) an on-axis read-out of a reflective phase mask (51) wxth 
a beam splitter (52) . Both configurations (A, B) may be 
utilized in the systems shown in Figs. 1-3. 
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PHASE ENCODING FOR DC PHASE FILTERING 

In the following an example of encoding a spatial phase 
maS K and a spatial phase filter will be given based on a 
system filtering in the DC-frequency range. The system 
chosen in this example is based on a 4-f lens 

in Pin i =nd illuminated by 
configuration as shown in Fig. 1 ana 

^- in the visible frequency 

electromagnetic radiation in tne 

domain, hereafter simply denoted as light radiation. 

Assuming that the illuminating light is monochromatic and 
has a substantially flat amplitude profile we obtain the 
following spatial amplitude distribution emitted from the 
spatial phase mask: 

*x.y) = rect(£,£)exp(i*x.y)) 



where «<*. y) = V» "P" sent the SPatiaUy 

phasor values and Ax4y is the area of the input phase 
20 modulating spatial light modulator. 

It turns out to be convenient to separate a(x, y) .into two 
terms describing a spatially invariant DC-value, a . and a 
spatially varying AC-contribution Aa(x, y) . The DC-value 
25 can be found as: 



(10) 

a = - ~ 
AxAy 



JJexp(i<Kx, y))dxdy 



30 



Subsequently the AC-term is expressed by: 

Acrfx, y, = «p(i*x. yfl - ^ fl.*"* >•>*<* « > 1 
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The separation of «* y) a spatially invariant DC- 

term and a spatially varying AC-term is an important 
; o^t and will he used throughout the remaining part o 
5 this example, especially in the descrrpUon of the 
spatial filtering procedure. 

The spatial filter utilized in this example is chosen as 

circular phase contrast filter .different transver 

,0 shapes can also he used, centered around orrgo rn the the 

„ j nma in denoted by coordinates [f x , f y ) • 
spatial frequency domain, oenoue * 



2tf r ) = 1 + (exp(itf) " * iT ijt) 



(12) 



„ where f, - Jtt ^ ^ 

V, describes the size of the circular (circ, phase 



filter. 



20 



25 



30 



In the spatiaf frequency domain (the filtering plane, the 
Pour er transformation ,3, of the spatially modulated 
Fourier t soatia l phase mask is present, 

licrht radiation from the spatial pim 

iignt i-av Frmrier transformed light 

The filtering operation on the Fourier 

Ration performed by ^ -££^1 

Utered li«ht is inverse Fourier transformed by he 

filtered ixg (Fourier transformation and reflected 

second Fourier lens (Fourie amp litude 
output coordinates) and the resulting p 

n i an p (with coordinates , y ; i 
distribution in the image plane (witn 

can accordingly be written as: 
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o(x', /) = /) * («M - y))cir fe)) 

= a(x', y') + a(exp(i^) - , ^) <* 3 > 

= [exp(i^x', y')) + «(exp(i^) - * ect (f^ ' Ay") 

fv' v'1 g 9T . outlined by 
Within the illumination-region, (x,y)e.K, 

rectf^^V ° ne obtain5: 

5 H ,, yi .x + fl.^gH.i<g-^.-*-'" l+ S] 1141 

t-h*t lofx' y'f = 0 corresponding to complete 
Requiring that |o\x„, y c ;| 

the lowest intensity level m 
darkness as tne 

,0 regions (x' 0 , y' a ) € «; implies: 

/A - ^x'.v'l has been used, 
where the abbreviation 0 e - ^x 0 , y 0 ; 

The solutions to Eq. (15) are given by: 



sir{ h - 4.*j\±^^±Jtl (16) 

R- — — - 



2 sin 



(!) 



20 The requ 



requirement o < R < 1 ^P iies that: 
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+ I** P = 0,±1,±2, . • . 



leading to 

^ i < ft < 1 ( 18) 

2 sin 



" in llJ 



5 whe re the + sign is for Rvalues in the interval: 



10 and the -sign is for Rvalues: 

The corresponding interval for (<f>- a - t.) ^ 



5*r . _ o» (1 9 ) 



15 



(#5 -#J e ]f ; 'f[ 



(20) 



(21) 



inserting the expression for R , one 
intensity expression: 

20 



obtains the simple 



|0(x , yi - si^ - y') * f)] 

where l-l -I j.\ (23) 

J| exp(i «x, y))dxdy = AxAyR exp(a * ; ) 
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The phase-only translations imply that eneroy is 
conserved: 

(24) 

3V*# y )\ - — ■> 



A special case: 

The most convenient choice for 5 is: a - - W*** 
tha t • - -PU-2-.. so that the output intensity can he 
10 described as: 

|o(x'. yf = 2[l - cos(^'. y'»] (25) 

xn this case the phase^intensity *appin 9 is describee by 
15 the intervals [0; a] -* [0*] ■ 

By settin, a - \ one obtains the followin, reouire.ent to 

the phase function #x, y) : 

AxAy 

JJ cos(^x, y))dxdy = {26) 

AxAy 

JJ sin(<*(x, y))dxdy = 0 



20 



• „ f«r Ux' v'tf in Eq. (24) yields: 
inserting the expression for |o(x . y J 

2 {J [1 - co^x', y'))]dx'dy' = AxAy <* 7 > 



AxAy 



., h the first of the integral expressions 
■>5 in accordance with the nrst 



♦5 

in Eq. (26) 
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- Encoding procedure: 

. A given intensity distribution (image) |o(x , y )| 
desired at the output side of the optical setup. 

5 . Pixellation of the image, that is generally 
5 Resented in the greyscale range: [Cgmax], P-v.des 

the relation: 
JJ|o(x'.yfdx'dy' = AxAy => EK*- * " ^ p1 ^" ' 



10 • The 



for the desired image |o(i, jf 1S 
histogram for tne 

adjusted (ad3) within the greyscale range [Cgmax], so 
that the previous point is fulfilled: 

|o(i, j)f -> |°U • 



„ . The phase values can now be calculated as: 



Hi, j) = arccosjl — 

. As before pixellation provides the relation: 

£ sin(#i, j)) = 0 . 

20 w filled by complex 

. The previous point can now be y >m 
conjugating half the input pixels having 
phase value in the phase histogram. 



25 



- nh.se flipping provides a valuable 
The phase ^"^//JL, for manipulating the 
tool (an extra degree of ^ [he 

spatial frequency content in 
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separation of low and high frequency terms at the 
filter plane. 

The scheme is robust to constant phase errors across 
the input spatial phase modulator, since Eq. (22) is a 
function of the difference: ^-(Kij). only. 
Furthermore, small variations in the individual pixel 
phase values do not introduce any detrimental effects 
because the average value a , is a result of a very 
large phasor sum. 



. If the desired intensity distribution is too small to 
include all energy, that is, the histogram is scaled 
to maximum and the left hand side of Eg. (24, is still 
smaller than the right hand side, then the input phase 
object can be scaled until Eg. (24) is fulfilled. In 
order to obtain a scale invariant output intensity 
level a dynamic focusing system is needed. Similarly, 
intensity invariance can be obtained by controlling 
20 the radiated power from the light source 

Alternately, one can ignore the residual background 
illumination and obtain intensity levels with a gain 
factor of 9" (background constant equal to 1 ) for 
narrow generally shaped line structures (e.g. Eq. 

25 (14)). 

Example 1. H the individual steps 

A very simple example illustrating tne 

in the above procedure will be given below. To simplify 
,0 the example it will be considered in one dimension only 
The starting point for encoding the spatial phase mask in 
this example is based on the following parameters: 



WO 96/34307 



28 



PCT/DK96/00190 



a = 0.5 

$ = n (28) 

grin ax = 4 

■„ fh. oixellated 3-step function shown in Fig. 5 

Consider plane as an inte nsity 

to be synthes.zed in ^ ^ 

5 :rrr -.s ~ ^ :r in 

spatial phase mask and the image intensity values. 



|o(i)f = 2[l - coa(tfi))] 



10 



To proceed 
accumulated 



from here 
intensity 



it necessary to 
2 \<*if in the 



(29) 

calculate the 
image to be 



15 



. , 9md The accumulated intensity is easily 
synthesized The histograro whe re the x-axis 

calculated from an image h g represents the 

arevlevel value and the y-axis re P i.« 
represents greyleve grey level value, 

amount of pixels in the image at a given g 
By US e of a histogram ? ^ simply found as the 

- all areY level values (x-axis) multiplied 
weighted sum of all greyle describ es, so to 

r^rca^ations are not needed since we only have 
3 greylevels with well-defined separations. 

The value for the accumulated intensity has to obey the 
25 equality: 



20 



gmax n 



(30) 



From Fig. 5 w « obtain: 
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£ |o(i)| 2 = 4pixels • 0 + 4pixels . (0.5max) + 6pixels • max = 8max 



(31) 



5 So 



that the value for max can be estimated to be: 



10 



15 



20 



7 (32) 
max = - v 

4 



The corresponding adjusted intensity levels, \c(i)\\ di , are 

therefore: 7/4, 7/8 and 0. These values can now be 

utilized to calculate the phase values of the spatial 
phase mask from the relation: 



25 





( mi,} 

= arccos 1 ■ ■ 


gmax 

V > 


I ) 



(33) 



where from we obtain the three phase values: 1.45 rad. 
0.97 rad. and 0 rad. 

The last step needed in order to encode the spatial phase 
mask is that the following equality is fulfilled: 



£sin(tfi)) = 0 



(34) 



Since we have the choice to use complex conjugate phasor 
values (two phasors giving the same intensity level) many 
approaches can be taken from here. A simple approach is 
to flip every second phasor with its complex conjugate 
value as shown in Fig. 6. The final phase values used in 
the phase mask are accordingly: ± 1.45 rad. ± 0.97 rad. 
and 0 rad. 



30 
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As the last step we can check whether the criteria: 
J. 1/2. is actually fulfilled with the chosen phasor 
encoding: 

14 (35) 



WO 96/34307 



31 



PCT/DK96/00190 



• GENERAL PHASE CORRECTION PROCEDURE INTEGRATED WITH THE 
PHASE ENCODING 

In Eq. (14) we obtained an analytic relation between the 
5 phase values in the spatial phase mask and the resulting 
intensity distribution, within the region (x', y') e SR' : 



0 The analysis leading to the above relation was based on 
the assumption that ft is a constant value within the «»- 
domain. In other words, the following approximation was 
applied: 

3-(s«* y)^{h-]\ s ~ ect (£ ' 1371 

IS 

However, for certain spatial filter parameters the 
lefthand side of this expression will not be a space 
invariant constant value throughout the whole *' -domain 
but will instead manifest slowly variations/oscillations. 
20 This will introduce small errors in the final 
superposition between the phase filtered DC-value and the 
direct propagated AC-signal. In order to circumvent this 
problem a technique is needed that can counteract the 
distortions by use of phase-only encoding in the 

25 components already present in the system. In what follows 
a procedure for integrating predistortion that 
counteracts the above mentioned distortions will be 
described that is purely based on modifying the phasor 
values in the spatial phase mask at the input side of the 

30 system. The method can also counteract other types of 
distortions inherent in a practical implementation of the 
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.h.rmore the method can be applied in systems 
system. Furthermore, cne 

filtering at other spatial frequences than DC. 

Procedure: . £ ^ help£ul to 

^ encoding the « P ^ ^ input phase 

h3Ve tion as a function of an adjusted (electronic, 
distribution as a addre ssing the input 

image grey- level distribution. I.„. addressing 

spatial light modulator: 

ginax L (38 ) 



«^ == 1 + 



k « taken into account that a(x', y') is not 
where it has been taken into o 

dered as a constant but manifests a smooth 
IS considered as a image domain> 

oscillating behaviour within the P 
The maximum value of I sln is denoted gmax . 

r*n derive a formula for the grey-level 
NOW ' In Tx (x< '), that one needs t o apply in order 

20 C ° rreCtl " a l se function that compensates for the 
to encode a phase ^ 
s patia! variation of the average phase value a(x , y ) = 



f ^.ix*. /)^:,,(x-,y'))| 

#x\ y') = acosl 1 - J 



(39) 



25 



wher e the second relation has been derived from the first 
by setting a - 1 / 2 and 6 = « • 
. B y .nserting the second Nation :n the first express., 
one qets: 

i 
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Ai sJn (x', y') = 



- 1 



[2|S(x', y')| J 



(40) 



5 This formula is however not directly useful because it is 
related to the histogram adjusted grey-leve! distribution 

denoted by I slm . 

One needs a formula that relates the above correction 
,0 term to the original input grey-level distribution l(x, y) 
that has not been modified by histogram adjustments. This 
ls important since the effect of the grey-leve! 
corrections also have to be incorporated in the procedure 
of histogram adjustments. 

15 

The histogram scaling gives: 



l(x, y) = y—- *s>»(x> y) 

max 



(41) 



20 where J m and are the maximum grey-level values 

occurring in the original and the adjusted electronic 
grey-level distributions respectively. 

Similarly, one can apply this relation to the intensity 
25 correction term AJ,„ and obtain: 

««. y) = I(x. y) * al(x. y) - -j>— (!„.<* y) + **...<* *» 

(42) 



resulting in: 

JO 
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A ' Y ' 2ja(x, y)| I v 

,„ order to have enough dynamic range In grey-levels for 
C he Erection term one can derive an ineguality from the 
, ab ove relation by using the fact that I.. * v™ax : 



3 



Cm V 5j "" nax 



(44) 



or 



2|a»in| • amax (45 ) 



10 



sl „ce the first term is the dominating term in the 

„, m for the intensity correction it will 

:n Efficient 3 ust to have the much simpler 

corrections : 



15 

l(x, y) 



-Tmax - 2 | a " ,l "l ' gm3X 
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PROPOSED APPLICATIONS: 

. La ser machining, marking, branding, trimming, hardening, 
5 scribe, labeling, welding and cutting on two and 
^ee-dimensional surfaces especially by use of C0 2 and 
Nd'YAG laser based systems. The main advantage is that 
7»Z is not absorbed in the system (thereby preventing 

«z of ^ ^-^T^j^rT^ 

delivered to selected regions on 
plane. High power can be delivereo to 

a work piece simultanously . 

I5 . cmcient and dynamic spot-array generators based on 
ohase contrast imaging. In order to provide bras or 
To d n, beams for arrays of electronic 
as bistable elements, photonic switches and smart prxels. 



20 



25 



* ^r-tnred light (lossless) for machine 

. ^"" OT i °\;" UCt E Ur g ed p r odic and sKew periodic mesh 
vision applications. E.g. pe" 

,rid illumination that can be updated in parallel. 

,. ^- ~ c /laser 3D direct writing 
. Photolithographic applications (laser 3D 

in parallel without the need for sequential scanning). 
in parallel ^ waveguides in Ge - 

E.g. high power laser direct 

doped silica. 

. spatia! light intensity modulation in genera! by use of 
pure phase modulation (radiation focusators, . 



. Laser beam shaping (dynamic I . 
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. Highly efficient parallel image projection without the 
need for a laser scanning device. 

. Dynamic Infrared Scene Projection (DIRSP) . 

. Exposure device for grating and mask production. 

• LIDAR applications. 

• Laserprinting in parallel. 

• Lasershow applications. 

• Atmosphere research. 
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CLAIMS 



1 A method for synthesizing an intensity pattern with low 
loss of electromagnetic energy, comprising spatial modulation 
of electromagnetic radiation with 
5 a spatial phase mask for modulation of the phase of the 

incident electromagnetic radiation by phasor values of indi- 
vidual resolution elements of the spatial phase mask, each 
phasor value being determined in such a way that 

1) the values of the Fourier transformed phasors attains 
10 predetermined values for predetermined spatial fre- 
quencies , and 

2) the phasor value of a specific resolution element of 
the spatial phase mask corresponds to a distinct 
intensity level of the image of the resolution 

15 element in the intensity pattern, and 

a spatial phase filter for phase shifting of a part of 
the electromagnetic radiation, in combination with 

an imaging system for generation of the intensity pattern 
by interference in the image plane of the imaging system 
20 between the part of the electromagnetic radiation that has 

been phase shifted by the phase filter and the remaining part 
of the electromagnetic radiation. 

2. A method according to claim 1, wherein the fact that for 
each resolution element of the spatial phase mask, each 
25 intensity level of the synthesized intensity pattern may be 
represented by one of two different phasors of complementary 
phase values is used to control the phase of the Fourier 
transform of the phasors at specific spatial frequencies by 
selection of phasors with appropriate phase values in such a 
30 way that the function defining the intensity values as a 

function of the phasor values comprises two regions of biuni- 
que functional dependence between phasor values and corre- 
sponding intensity values. 
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i i~ i rs-r o wherein the modulus of 
3. A method according to claim 1 or 2, wherein 

the Fourier transform o£ the phasors at specific spatial 

regencies is adjusted in order to control the range of 

intensity levels of the synthesized intensity pattern. 

5 4 A method according to any of claims 1-3, wherein the phase 
5 of Courier transform of the phasors at specific spatial 
frequencies is adjusted in order to control whether the 
relation between each phasor and the corresponding intensity 
level is a monotonic increasing or a monotonic decreasing 
10 function. 

5. A method according to any of claims 1-4, 
adj ustment of the modulus of the Fourier transform of the 
phasors at specific spatial frequencies comprises at least 
one of the following measures: 

». adjusting the individual phasors of the resolution 
lements of the phase mask maintaining prescribed relative 
elements of the pna . n ^ naities of resolution elements 

intensity levels between intensities ot res 

of the intensity pattern, 

b , adjusting the individual phasors of the resolution 
20 elements of the phase mask by histogram techniques, 

spatially scaling the phasor pattern of the phase 



c) 
mask, 



25 



d) utilizing half tone coding techniques. 

6 a method according to any of claims 1-5, wherein 

of L radiation source is controlled in response to the 

intensity range of the intensity pattern. 
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,. A method according to any of claims """f^ 
„hasor of the phase mask is selected from the set of two 
STors with elementary phase values in such a ^y hat a 
specific spatial frequency distribute of the ntensity 
5 L electromagnetic radiation at the spatial filter 
attained. 

8 a »ethod according to claim 7. wherein the phase of phase- 
's of adjacent resolution elements and of substantially 
identical phasor values alternates between the two possible 

10 complementary phase values of phasors. 

9 A method according to claim 8. wherein the shape of the 
pbaseTilter is adapted to match the spatial frequency con- 
tent of the phasors of the spatial phase mask. 

10. A method according to any of claims 1-9. >^££Z. 
15 of the phase mask and the phase shift of the phase filter are 
selected in such a way that the intensity level of at least 
"e resolution element of the synthesized intensity pattern 
is zero. 

U. A method according to any of claims 1-10, wherein the 

» — - «- f °::Ts ssr s rirro^rte 

TZ\o obtain a maximum dynamic range of the phasors of 
the phase mask. 

25 12 A method according to any of claims 1-11. wherein the 
pLe filter is centred around the -reorder diffraction 

region. 

^n^-im 19 wherein at least one 
13 A method according to claim 12, wherein « 

intensity level of the intensity pattern may be represented 
30 by one of two complex conjugate phasors. 
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M A method according to claim 12, wherein at least one 
density level of the intensity pattern is represen t e d hy 
one of two phasors with phases the sum of whrch equals 

U A method according to any of claims 10-14. wherein the 
ILL of average value of the phasors and 
introduced by the phase filter of the 

electromagnetic radiation are related to each other u> accor 
dance with the equations 



HI--T-TT7 and |a|Sl 

2 | sin- | 

10 16. A method according to claim 15. wherein modulus of the 
average value of the phasors ranges from 0.1 to 0.9. 

„. A method according to claim 15 or 16. wherein modulus of 
the average value of the phasors ranges from 0.25 to 0.75. 

1.. A method according to any of claims 15-17. wherein 
„ modulus of the average value of the phasors ranges from 0.4 
to 0.6. 

19 A method according to any of claims 15-18, wherein 
^Jof the average value of the phasors is appro X1 mately 



0.5 



2 „ 20 . A method according to any of claims 1-19 «^ein the 
pnase shift of the phase filter ranges from ,/4 to 7r/4. 

21 . A method according to any of claims 1-20. wherein the 
phase shift of the phase filter ranges from ,12 to 3r/2. 

22. A method according to any of claims 1-21 wherein the 
25 phase shift of the phase filter ranges from 3»/4 to 5./*. 



10 
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23 A method according to any of claims 1-22, wherein the 
phase shift of the phase filter is approximately ». 

24. A method according to any of claims 1-23, wherein the 
phase shift of the phase filter is variable. 

5 25. A method according to any of claims 1-24, wherein the 
imaging system comprises zooming means for variable scaling 
of the intensity pattern. 

26. A method according to claim 25, wherein the zooming of 
the imaging system is dynamically controllable. 

27 A method according to claim 25 or 26, wherein the zooming 
of the imaging system is controllable in dependence of the 
scaling of the phase mask. 

28 A method according to any of claims 25-28, wherein the 
power of the radiation source is controllable in response to 
the spatial scaling of the pattern in the phase mask and/or 
the zooming of the focusing system. 

29. A method according to any of claims 1-28, wherein the 
phase mask is a spatial light modulator. 

30. A method according to any of claims 1-29, wherein the 
optical function of a Fourier- transforming lens is encoded 
into the phasors of the phase mask. 

31. A method according to any of claims 1-30, wherein the 
optical function of an output lens is encoded into the phase 
filter. 

32. A method according to any of claims 1-31. wherein part of 
the electromagnetic radiation modulated by the phase mask has 
a substantially flat intensity profile in the image plane. 



15 



20 



25 
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■ ,3 a method according to any of claims 1-32. wherein the 
JTt electromagnetic radiation comprises one or more 
U U gh sources o di «erent wavelengths corresponding to three 
a Mere" colours, such as red. green and blue, for gener- 

5 ation of intensity pattern of arbitrary colours. 

34 A phase contrast imaging system for synthesizing an 
intensity pattern with low loss of electromagnets energy, 
comprising 

at least one source of electromagnetic radiation for emission 
10 of electromagnetic radiation. 

a spatial phase masK for modulation of the phase of the 

eachVasor value has been determined in such a way that 

the values of the Fourier transformed phasors attains 
predetermined values for predetermined spatial fre- 
guencies, and 

2) the Phaser value of a specific resolution element of 
the spatial phase masK corresponds to a distinct 
intensity level of the image of the resolution 
element in the intensity pattern, and 

a spatial phase filter for phase shifting of a part of the 
electromagnetic radiation, and 

f „r- operation of the intensity pattern by 
an imaging ^° r J™ o£ the iina9i „g system between 
interference in the image plane or t * 
the part of the electromagnetic radiation that has be P 
shifted by the Phase filter and the remaining part of the 
electromagnetic radiation. 



15 1) 
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• 35 A system according to claim 34, wherein the imaging 
system comprises a first and a second Fourier transferring 
lens, the spatial phase mask being positioned in the front 
focal plane of the first lens, the spatial phase filter bexng 

5 positioned at the back focal plane of the first lens, and the 
second lens being positioned so that its front focal plane xs 
positioned at the position of the back focal plane of the 
first lens. 

36 A system according to claim 34 or 35, wherein the imaging 
10 system comprises one Fourier transforming lens the spatial 
phase filter being positioned at the back focal plane of the 



lens. 
37 



A system according to any of claims 34-36, wherein the 
imaging system comprises one imaging lens the spatial phase 
15 filter being positioned in the back focal plane of the lens. 

38 A system according to claim 34-37, further comprising a 
polarizing beam splitter and a quarter wave plate and/or a 
phase filter reflecting electromagnetic radiation incident 
upon it. 

20 39. A system according to any of claims 34-38. wherein the 
spatial phase filter changes the phase of the DC-part of the 
electromagnetic waves and leaves the phase of the remaining 
part of the radiation unchanged. 

40 A system according to any of claims 34-38. wherein the 
25 spatial phase filter do not change the phase of the DC-part 
of the electromagnetic waves and changes the phase of the 
remaining part of the radiation. 

,1 A system according to any of claims 34-36. wherein the 
spatial phase filter blocks the DC-part of the electromag- 
30 netic waves and leaves the refining part of the redrawn 
unchanged . 
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« A system according to any of claims 34-41. wherein the 
source of electromagnetic radiation is a Laser. 
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